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I. INTRODUCTION 

Considering theories with massive gravitons |lj is motivated by the observation of the 
current acceleration of our universe |2j, since the graviton mass can induce an effective 
cosmological term. Although such theories typically exhibit unphysical features, as for 
example the Boulware-Deser ghost J^j], the recent discovery of the special massive gravity J4| 
and its bigravity generalization Q which are ghost-free p] suggests that such theories can 
indeed be good candidates for interpreting the observational data. This motivates studying 
cosmological solutions with massive gravitons. 

The first self-accelerating cosmologies in the ghost-free massive gravity were obtained 
without matter, they describe the pure de Sitter universe l7fl. The matter source was then 

n 

included for special values of the theory parameters |8|. A surprising feature of these solu- 
tions is that their physical and fiducial metrics do not share the same Killing symmetries - 
although the physical metric is of the Freedman-Robertson- Walker (FRW) type, the fidu- 
cial metric is inhomogeneous. This creates a technical difficulty, since the equations for the 
Stiickelberg fields reduce to a non-linear PDE. 

One more similar solution was found in Ref.jsl, where it was argued that, even though 
the background geometry is homogeneous and isotropic, the inhomogeneous structure of 
the Stiickelberg fields should be visible when the background is perturbed. However, this 
effect should be strongly suppressed by the smallness of the graviton mass, so that in small 
compared to the graviton Compton wavelength regions the deviations from the standard 



FRW cosmology should be small. T 
cosmologies, since, as was noticed in 



'hese properties seem to be generic for massive gravity 
the 'genuinely' homogeneous and isotropic solutions 
for which both metrics would be FRW do not exist in the theory, at least when the metrics 
are spatially flat. Even though such solutions were later found in the spatially open case, 
they are less interesting physically, and show in addition a nonlinear instability 10j. 

It seems therefore that physical cosmologies with massive gravitons should be described 
by solutions with a homogeneous and isotropic physical metric, but with an inhomogeneous 
or anisotropic (or both) fiducial metric (unless the latter is chosen to be non-flat Jll|). 



Such solutions were constructed also in the bigravity context |12|, when both metrics are 
dynamical but are not simultaneously diagonal and do not share the same symmetries. 
However, an all cases the solutions were obtained only for constrained and not generic 
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values of the theory parameters. Very recently, a solution for generic parameter values was 



announced in 



15| . but without determining the Stiickelberg fields, so that the most difficult 



)art of the problem was actually skipped. Finally, the complete solution was obtained in 
16| . both within the bigravity and massive gravity. In what follows we shall present a 
detailed derivation of this result. 

For pedagogical reasons, we shall restrict our discussion below only to the massive gravity 
case, since in the bigravity the procedure is essentially the same but the formulas are more 
complicated. The analysis of (l6j was carried out within the tetrad formalism, but we shall 
employ below the standard parametrization in terms of the matrix square root used by most 
authors. At the same time, we shall show how the same results are recovered within the 
tetrad formulation. We shall try to be maximally explicit. For example, it turns out that 
it is not easy to find in the literature the explicit form of the field equations in the theory, 
and most authors prefer to put a symmetry ansatz to the action and then vary. We shall 
therefore show how to vary the action in the general case. We shall also show that the tetrad 
formulation is equivalent to the standard one. 

We make the ansatz for which the physical metric is FRW but the fiducial metric is only 
spherically symmetric and not diagonal. We then calculate the matrix square root and show 
that the resulting equations are such that the Stiickelberg scalars effectively decouple, and 
the metric satisfies Einstein equations with an effective cosmological term plus the matter 
source. However, such a decoupling is only possible if the scalars satisfy a consistency 
condition expressed by a non-linear PDE. Fortunately, the latter can be solved exactly. As 
a result, we obtain the most general cosmological solution for which the physical metric 
is homogeneous and isotropic but the Stiickelberg fields are inhomogeneous. The solution 
includes the matter source, it exists for generic parameter values, and it describes a FRW 
universe that can be spatially flat, open or closed, and which shows the late-time acceleration 
due to the effective cosmological term mimicked by the graviton mass. 

The rest of the paper is organized as follows. In Sections II, III we introduce the ghost-free 
massive gravity and derive its equations of motion. The symmetry reduction is described 
in Section IV. The solution for the metric is presented in Section V, while the consistency 
condition for its existence is analyzed in Section VI. In Section VII all results are rederived 
again within the tetrad formulation. Finally, in the Appendix we describe all possible 
solutions for which both metrics are homogeneous and isotropic. 
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II. THE GHOST-FREE MASSIVE GRAVITY 

The ghost-free massive gravity theory of de Rham, Gabadadze, and Tolley (dRGT) Q is 
defined on a four-dimensional spacetime manifold equipped with the metric g^ u and carrying 
four scalar fields X A (Stuckelberg scalars) which parameterize the fiducial metric f^ v = 
VABd fl X A d u X B with t]ab = diagfl, —1, —1, —1]. The action is 



(2.1) 



where m is the graviton mass and S m describes the ordinary matter (as for example perfect 
fluid) that interacts with g^ v in the usual way. To define U, the interaction between g^ v and 
f^v, the key element is the tensor 



7t = Vd^fau, 

where g^ v is the inverse of g^ v . Here the square root is understood in the sense that 

(7 2 r„ = Y a ll = g^fav, (2.2) 

or, using the hat to denote matrices, 

f = g-'l (2.3) 

Introducing = 6% - i% with traces [JC] = tr(/C) = and [/C n ] = tr(/C n ) = (/C n )£, the 
interaction is 

U = U 2 + a 3 U 3 + a A U A , (2.4) 

where a 3 , a A are parameters, and 

1 

2! 
1 

3>. 
1 



W 2 = -([/C] 2 -[K 2 ]), 
W 3 = ^([/C] 3 -3[/C][/C 2 ] + [/C 3 ]), 



U A = - ([/C] 4 - 6[/C] 2 [/C 2 ] + 8[/C 3 ][/C] - 6[/C 4 ]). (2.5) 
Equivalently, with Xa being eigenvalues of /C, 

U 2 = X a X b, Us = X ^b\c, U a = X^X-sK = det(/C), (2.6) 

A<B A<B<C 
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One more equivalent representation is 17 1 



^2 = - y e^ upa e at3pa JC£ICp, 

u 3 — 7y tfivpcrt /V a /V -/3' V 7) 



u, = -1 . (2.7) 



[one has e 0123 = — e i23 = 1). 



III. FIELD EQUATIONS 

Some care is needed when varying the action with respect to the metric. Let us first vary 
the constraint ( 12.2ft . 

Sjj + j6j = 6g~ x f, (3.1) 

which gives 

.jHW^^/r 1 . (3.2) 

This cannot be resolved with respect to Sj, because 7 and 5^ do not commute. However, 
only the matrix traces enter the action. Taking the trace of (13. 2p gives 

tr(5 7 ) + tr^yy -1 ) = 2tr(5 7 ) = tr^W 1 ), (3.3) 
and noting that / = gg~ l f = gj 2 , one obtains 

6[K\ = -tr(5^) = ~tr(5g- 1 gy). (3-4) 

Next, one has 

5[K, n+1 ] = 5tr(/C™ +1 ) = (n + l)tr(<5/C/C n ) = -(n + l)tr(5 7 /C n ), (3.5) 



where the cyclic property of the trace is used. Now, multiplying ( 13.2ft by K7 1 from the right, 
taking the trace and using the fact that 7 and AC commute, one obtains 

2tr(5 7 /C")=tr(5r 1 /V 1 /C"), (3.6) 

which finally gives 

8VC+ 1 ] = -^triSg- 1 ^) = VWU£"K- (3.7) 
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It turns out that the matrix g^KJ 1 is actually symmetric 13|] . To see this, one uses the 
relation (m is a non-negative integer) 



gNt 1 /) g-'-^ft 1 ) , (3.8) 

which follows from the fact that squaring it gives identity. This implies that 



9 ( Vr 1 / ) =[\/f9- 1 ) 

and hence 



m \ tr / / , \ m \ tr 

tr 



J (3.10) 

so that g-j m is a symmetric matrix. Since g^K n can be represented as the sum of terms of 
the form g^ m with different m, it is also symmetric. This finally gives 

W n+1 ] n + 1 - 

~s^r = — — )m > ( 3 - n ) 

and the expression on the right here is symmetric with respect to \i •f-)- v. 

It is now straightforward to vary the action. This gives the Einstein equations 

G^ = m 2 T tlu + MGT^\ (3.12) 

where Tfu, is the matter energy- momentum tensor obtained by varying S m , while 

T»» = 2-^L-Ug» v . (3.13) 

Using U in (12. 4p , (I2.5P and applying (13.1 ip gives 

-a 3lfia {U 2 5 a u -\K]K, a u + (K 2 ) a u } 

- a 4l , a {U 3 5 a u -U 2 K a u + [/C](/C 2 )? - (/C 3 )?} 

-Ug^, (3.14) 

where ^ m = g^ a l a a - The above considerations guarantee that 7) tJy = T^ v y It is also worth 
mentioning the equivalent representation, 

T,* = lim {K - [K]8« + a 3 Z« + a^} -U 9fiu , (3.15) 
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where 

^ = ^ e^e^/C^/q. (3.16) 

Since the matter energy-momentum tensor is conserved due to the diffeomorphism-invariance 
of the matter action S m , 

the Bianchi identities for the Einstein equations imply the conservation condition, 

= 0, (3.18) 
which can be viewed as equations for the Stiickelberg fields. 

IV. SYMMETRY REDUCTION AND TAKING THE SQUARE ROOT 

Let us now choose spherical coordinates x^ 1 = (t, r, tp) and assume the physical metric 
to be homogeneous and isotropic, 

ds] = &(t) 2 dt 2 - s?{t)dr 2 - R 2 (t,r)(d$ 2 + sin 2 ??V), (4-1) 

with R(t,r) = a(t)/k(^)- Here /k = {r, sin(r), sinh(r)} for k = 0, 1, —1, which corresponds, 
respectively, to spatially flat, closed, or open FRW universe. The metric is invariant under 
spatial rotations and translations. We also assume the matter to be a homogeneous and 
isotropic perfect fluid, so that 

8nGT {m)p x = diag[p(t), -P(t), -Pit), -Pit)]. (4.2) 

The conservation condition ( 13. 1 7[) then reduces to 

p + 3-(p + P) = 0. (4.3) 
a 

As for the fiducial metric, we assume it to be only spherically symmetric, but not necessarily 
homogeneous, 

ds) = dT 2 {t,r) - dU 2 {t,r) - U 2 {t, r){d$ 2 + sin 2 M^ 2 ) = (4.4) 
= (T 2 - U 2 )dt 2 + 2(TT' - UU')dtdr + (T /2 - U l2 )dr 2 - U 2 (t, r){d$ 2 + sin 2 #d<p 2 ). 



s 



Therefore, the two metrics do not have the same Killing symmetries. 
The above expressions imply that 

( A C \ 
-C B 



(7 2 r, = <r/< 



U 2 /R 2 
\ U 2 /R 2 J 



where 



A _ t 2 - U 2 TV - UU' 



U 12 - V 2 



B 



1", 



To take the square root of this matrix, one makes the ansatz 

(a c 0\ 
-c b 
!i 
\ u ) 

Inserting this to (14. 5 j) gives algebraic equations 

a 2 - c 2 = A, 
c(a + b) = C, 
b 2 -c 2 = B, 
u 2 = U 2 /R 2 , 



whose solution is 



with 



A + A 
~^Y~' 



B + A 



C 



u 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



Y = VA + B + 2A, A = vMfi + C 2 . 



(4.10) 



V. SOLUTION FOR THE METRIC 

Since 7^ is known, we can compute the energy-momentum tensor in (|3.14j) . It is con- 
venient to lift one index and consider = g^ a T au . To begin with, the eigenvalues of 
K£ = 5£- Yu are 



2 - a - b ± J (a - b) 2 - 4c 2 , 
Ai, 2 = ^ , A 3 = A 4 = m, (5.1) 



which can be complex-valued. However, the symmetric polynomials in (12. 6p are always real, 

U 2 = u{u + 2a + 26 - 6) + c 2 + ab - 3a - 36 + 6 , 

U 3 = a3 (i - u )l( a + b- 2)u + 2c 2 + 2ab - 3a - 3b + 4] , 

W 4 = a 4 (l - «) 2 (c 2 + ab-a-b+l). (5.2) 

Inserting this and (14.71) into (I3.14p gives the following non-zero components: 

T ° = u (6 - 2b - u) + 36 - 6 + a 3 (u - l)(6u - 2u - 36 + 4) + a 4 (u - 1) 2 (6 - 1), 
T r r = u(6 — 2a — w) + 3a — 6 + a 3 (u — i)(au — 2u — 3a + 4) + a 4 (u — l) 2 (a — 1), 
T r ° = — c[(a 3 + a 4 )w 2 — 2(a 4 + 2a 3 + l)u + 3 + 3a 3 + a 4 ] , (5.3) 
T$ = T£ = u(3 — a — 6) — 6 — a6 — c 2 + 3a + 36 + a 4 (w — l)(c 2 + a6 — a — 6 + 1) 
+ a 3 [u{c 2 + ab + 3 - 2a - 26) - 4 - 2a6 + 36 - 2c 2 + 3a]. 

It turns out that 

c(T ° - T r r ) = (a - 6)T r ° . (5.4) 

The field equations G^ = m 2 Tjf + 8irGT^u imply that Tjf should be diagonal, since both 
G% and T^t are diagonal. Therefore, one should have 

T r ° = 0, (5.5) 

and so Tq = T r r if c ^ 0. Now, T r ° in (15. 3p will vanish if either c = 0, or if the expression 
between the brackets vanishes. The c = case will be considered in the Appendix. If c 7^ 0, 
then T r ° in ( 15. 3p will vanish if 

P 2 (m) = («3 + « 4 )^ 2 — 2(a 4 + 2a 3 + \)u + 3 + 3a 3 + a 4 = 0, (5.6) 

which requires that u is constant, 

_ 1 + 2a 3 + a 4 ± a/1 + a\ + a 3 - aj 
u— ■ , \p-() 

a 3 + a 4 

which will be real if a 4 < a 3 + a 3 + l. Inserting this to (15. 3p we find that the two components 

T ° =T r r = (1 - u)(u + ua 3 - a 3 - 3) (5.8) 
are also constant. The conditions V M T^ = reduce in this case to 

VX = 2 ^ W - It) = 0, (5.9) 
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which requires that Tq = T$. Now, using the above formulas we find that 

T ° - Tt = " + M <*3-2-a 3 [{u -a)(u-b) + c% (5.10) 

L — U 

and for this to be zero either the first or the second factor on the right must vanish. Let us 
require (the other option will be discussed in the next Section) that 

(u-a)(u-b) + c 2 = 0. (5.11) 

This condition guarantees that Tj^ is conserved, so that this is the equation for the Stiick- 
elberg fields. Assuming that this condition is fulfilled, T£ becomes proportional to the unit 
tensor and the field equations reduce to 

G p x = A5 p x + 8nGT {m)p x (5.12) 

with 

A = m 2 (l — u)(u + uaz — — 3). (5.13) 

The functions a, b, c effectively drop out and only the constant u remains. As a result, the 

effect of the graviton mass is the same as that of a constant cosmological term. Einstein 

equations (15. 121) further reduce to the Friedman equation 

sl ~\~ kci , . 

3 — = A + p, (5.14) 

Si 

where p(a) is determined by the conservation condition ( 14.31) . This equation describes a 
universe filled with matter and containing a cosmological term mimicked by the graviton 
mass. At early times the matter density p dominates, but in the long run the cosmological 
term wins, leading to the self-acceleration. 



VI. SOLUTION FOR THE STUCKELBERG FIELDS 

The above solution for the metric is only possible if equation (15. lip is satisfied, so that 
the whole procedure is consistent if only this equation can be fulfilled. Rewriting it as 

ab + c 2 + u 2 = (a + b)u (6.1) 

and using expressions (14. 9p for a, b, c yields 

(A + A)(B + A) + C 2 , 

— h u =uY. (6.2) 
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Now, using (14.1 Op one has 

(A + A)(B + A) + C 2 = AB + C 2 + {A + B)A + A 2 

= (A + £) A + 2A 2 = A(A + B + 2A) = AY 2 , (6.3) 

and so (16. 2p becomes 

A + u 2 = uY. (6.4) 

Next, using (14. 6p . one finds 



AB + C 2 = ; 

and therefore 



a 1 



(6.5) 



, TTJ' — UT' 

A = VAB + C 2 = , (6.6) 



while 



Y = VA + B + 2A = - V ( T + U ') 2 ~ ( T ' + ^) 2 ' ( 6 - 7 ) 
so that ( 16. 4p reduces to 



TV' - UT' + w 2 a 2 = way (T + [/' + T' + ?7)(t + U' - T' - U). (6.8) 

Squaring this finally gives 

{TU' - UT' + w 2 a 2 ) 2 - w 2 a 2 (T + U' + T' + U)(f + U' - T' - U) = 0. (6.9) 

This is a quadratic PDE for T(t,r), with the coefficients a(t) and U = uR = ua(t)f^(r) 
determined by the solution of the Einstein equation (15.141) . Although this equation looks 
complicated, some of its solutions can be obtained. 

Let us first consider the case of spatially flat universe, when U = ua(t)r. One makes the 
ansatz, 

T(t,r) = /(t)+Ca(t)r 2 , (6.10) 
with constant C. In this case Eq. (I6.9j) reduces to 

w 2 a 2 - 4Ca/ + 4C 2 a 2 = 0, (6.11) 

which can be resolved with respect to f(t) to give 

f a 2 . (u 2 



T(t,r)=C I -dt+l— + Cr 2 ) a . (6.12) 
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This solution agrees with the one obtained in Q for 03 = 04 = 0, when u = 3/2. 
When the universe is spatially closed and U = ita(t) sin(r), one assumes 



which reduces Eq. (16.91) to 



hence 



T(t, r) = f(t) + Ca(t) cos(r), (6.13) 



(C 2 + <)(a 2 + a 2 ) = / 2 , (6.14) 



T{t,r) = ±J \I{C 2 + w 2 )(a 2 + a 2 ) dt + Cacos(r). (6.15) 
When the universe is open and U = ita(t) sinh(r), one makes the ansatz 



reducing the problem to 



so that 



T(t, r) = f(t) + Ca(t) cosh(r), (6.16) 



(C 2 -n 2 )(a 2 -a 2 ) = / 2 , (6.17) 



T(t,r) = J V(C 2 - w 2 )(a 2 - a 2 ) dt + Cacosh(r). (6.18) 

This completes our construction, since we have determined the metric and the Stiickelberg 
fields for all three types of the universe and for generic a 3 , a 4 . 

The above solution exists if < 03+0:3 + 1. Let us consider the limit where this inequality 
is saturated, 04 = o 2 + 03 + 1. This implies that u + UO3 — 2 — 03 = 0, in which case the 
conservation condition ( 15. 9 j) will be fulfilled without imposing the constraint 05. 11 j) . This 
possibility has been first analyzed in [8| (and recently rediscovered in 14J). The obtained 
above general solution applies in this case too, and the physical metric is determined by 
the same equations (I5.14p as before. However, the function T(t, r) needs not to be now the 
same as before, since the constraint 05. 1 1 j) is no longer imposed, so that there is actually no 
condition for T(t,r). Therefore, T(t,r) remains arbitrary (a particular choice T = — j Udr 



was made in js|). 



VII. SOLUTION IN THE TETRAD FORMULATION 



Let us now see how the same results are obtained within the tetrad formulation used in 
16|. This formulation was originally introduced in 18| and then further developed in |l9| . 
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12|. 



it is equivalent to the standard formulation but can sometimes be more efficient 
Although originally it was used to argue in favor of a possible presence of the ghost |l8|, the 
absence of ghost in the tetrad description was shown in jl9|. 

The basic variables in the tetrad formulation are two tetrads e A and which determine 
the two metrics, 

g ^ = 7] AB e » A f» v = VAB ufa*, (7.1) 

and also the tensor 

ft = ■ (7-2) 
Defining K,^ = 5^; — 7^, the action is still given by Eq. fl2.ip . with the interaction 

U = \ ((/CP 2 - K%K$ - f e^TO - ^ e^e^K^K^ (7.3) 

(the notation C3 = —0:3, C4 = — was used in jisjl). Comparing with (l2.4p .( T2~Tl) . this is the 
same expression as before, up to the replacement 7^ — > 7^. Varying the action with respect 
to e A is straightforward and gives the Einstein equations = m 2 T^ + 8irGT^u with 

TZ = e^-U5». (7.4) 

Explicitly, 

Tpu = %a {K - [K-K + « 3 S: + a 4 Q a u } -Ug^, (7.5) 

where 7^ = g^vj^a an d an d defined by (I3.16p . This expression agrees with that in 
(I3.15p . up to the replacement 7^ — \ 7^. Therefore, if we could show that 7^ = 7^, this 
would mean that we have the same equations as before. 

The equality 7^ = 7^ can be enforced by the the following condition: 

UApf^B = > ( 7 - 6 ) 

with uj Ct i = Vca^, because 

T a Tu = efcd£eSa£ = e A ^ Aa e a B ^ = e A ^ Ba e a A ou* = g^f av , (7.7) 

and therefore 7^ fulfills the very same equation which defines 7^. The condition (17. 6p was 
originally postulated in 18|. However, it turns out it actually follows from the field equations 



fsee also 



19(|) . Indeed, the Einstein equations imply that T^ v is symmetric. The expression 



in (17. 5p will be always symmetric if only the first four terms on the right are separately 
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symmetric. Therefore, 7^ should be symmetric, and this guarantees that the other three 
terms are symmetric as well. As a result, the field equations require that 

9 m T v = 9vcn a ii • (7-8) 

Using the definition of 7^ and also g^ u = t\ab^ A ^, w here e A is the inverse of e^, these 
relations assume the form 

e J^cv = e°u c „ , (7.9) 

multiplying which by e A e B gives precisely Eq. (17.6jl . Therefore, the equality 7^ = 7^ is 

imposed dynamically, hence we can remove the tilde sign and conclude that the tetrad 

formulation gives the same theory as before. 

To obtain the solution, one makes the ansatz, 

Id Id Id Id 

a ot aar K dv Rsmv dip 



= a(adt + cdr), u 1 = a {b dr — c dt) , uj 2 = uRd'd ) u 3 = uRsm : ddip 



which fulfills (17.61) . with a = a(£), and where R = a(t)/k(r) is the same as in Eq.f l4.ip . This 
implies that g^ u is the same as in (14. ip . Computing 7 M V = e A u A then gives the same result 
as in (14. 7p . and therefore all analysis of Section |V] goes through without any changes. 

There remains to analyze the consistency condition (15. lip . Let us remember that the 
metric f^ y = T]AB^a u v should be flat. However, this does not mean that 1-forms u A = 
u^dx 11 coincide with differentials of the Stuckelberg scalars, because it is still possible to 
perform local Lorentz rotations, so that u A = L A B dX B where L A B is a position-dependent 
50(1, 3) matrix. Comparing with in (14. 4p . it follows that u 2 = Udd and u 3 = U sinfidip, 
while u° and u 1 can differ from dT and dU by a local Lorentz boost, 

tu° = cosh(a)c?T + sinh(a)<i£7, uj 1 = cosh(a)dU + sinh(a)dT , (7-11) 

where a is the boost parameter. Explicitly, 

a (a dt + cdr) = cosh(a)(Tdt + T'dr) + sinh(a)(£7 dt + U'dr), 

&(bdr -cdt) = cosh{a){Udt + U'dr) + sinh{a){fdt + T'dr). (7.12) 

Comparing the coefficients in front of dt, dr gives four conditions, which determine a, b, c, a in 
terms of T, U. Inserting the resulting a, b, c to (15. lip then gives precisely the same equation 
for T(t, r) as in (|6.8|) . Therefore, we recover the same results as before. 
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In the case when a 4 = a\ + a 3 + 1 discussed at the end of Section |VT| when the constraint 
(15. lip is not imposed, there is no condition for coefficients a, b, c, a obtained from (I7.12p . so 
that the choice of T(t, r) remains arbitrary. Equivalently, one can choose some value of a 
and then calculate T(t, r) from ( 17. 12ft . For example, setting a = requires that T' + U = 
and so T = — J Udr |8|] . 

VIII. CONCLUSION 

We presented the most general cosmological solution in the ghost-free massive gravity 
with flat reference metric for which the physical metric is homogeneous and isotropic but 
the Stiickelberg fields are inhomogeneous. The solution includes the matter source, it exists 
for generic values of the theory parameters, and its physical metric describes a FRW universe 
that can be spatially flat, open or closed, and which shows the late-time acceleration due to 
the effective cosmological term mimicked by the graviton mass. 

Even though the physical metric is homogeneous and isotropic, its perturbations are 



expected to be inhomogeneous, due to the Stiickelberg fields [9|. This effect will be propor- 
tional to m 2 and so will be small in small enough regions of space. However, it would be 
still interesting to compute the linear perturbation spectrum (see also j2Q|). It would also 
be interesting to see if the above construction could be generalized to describe non-linear 
anisotropic deformations of the homogeneous and isotropic background. 

We have also shown the equivalence between the standard parametrization of the dRGT 
theory and the tetrad approach (see also Jl9|). 

APPENDIX. SOLUTIONS WITH FRW METRICS. 

For the sake of completeness, we review in this Appendix solutions for which both the 
physical and fiducial metrics are homogeneous and isotropic. Most of them have been 
previously reported in the literature. They are obtained by setting c = in the formulas 
of Sections HVlVl Such solutions exist only for particular values of the spatial curvature k 
and/or show a degenerate fiducial metric. Therefore, they are less interesting physically 
than solutions described in the main text above. 

If c = 0, then the condition T r ° = in (15. 5p will be fulfilled without imposing the 
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constraint (15. 6p . while Eqs. (l4.8p will imply that C = 0, so that the metric f^ u in (14. 4p is 
diagonal. Eqs. (l^ . (0]l .( 0]l .( OU]l then reduce to 

f 2 -U 2 



2 _ 

U' 2 - V 2 



a 2 = - , (A.la) 



b 2 = -— , (A.lb) 



a 

TT' = UU', (A.lc) 

with U = na(t)/k(r). Here a,b,u are functions of t, r. The energy-momentum tensor is 
defined by Eqs. (15.31) with c = and should satisfy the isotropy condition, 

K = Tl (A.2) 

the homogeneity condition, (Tjf) 1 = 0, and the conservation condition, 

T ° + 3 - (T ° - T r r ) = 0. (A.3) 



The Einstein equations for the metric g^ v then reduce to 

3 = m 2 T ° + p. (A.4) 



Equations (IA. 1[) ( TA74|) determine a,b,u, a(t), and T(t,r). 



Let us consider first the isotropy condition (1A.2|) . Using Eqs. (l5.3p gives 

Tl — T$={u — b)[(au — 2a — 2u + 3)«3 + (1 — u — a + Gra)a 4 + 3 — a — u], (A. 5) 

and since this has to be zero, either the first or second factor on the right should vanish, 
which we shall call case I and case II, respectively. In the case I one has u = b and the 
conservation condition (IA.3[) reduces to 



- ka)P 2 {u) = 0, (A.6) 

with = u& and P2(u) defined in (15.61) . Depending on which of the two factors on the left 
vanishes, there are two subcases to analyze, let us call them la and lb. 
Case la: — ka = 0. One has a — /3/a and b — u — /3/a. Inserting this to (IA.la|) and 
(lA.lbp gives 



T' = Htf-l, T = pJ* 2 /h 2 + f 2 . (A.7) 
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These should fulfill flA.lcfl and also the integrability conditions df r T = d^ t T, which is only 
possible if — 0. Therefore, the Stiickelberg scalars are 



T(r) =pj dry/f*-l, 



U{r) = f3f k 



(A.8) 



with constant 0. The scale factor satisfies ( 1A.4I) with 



rriO 



„ 3/3(3a 3 + a 4 + 3) 3/5 2 (l + a 4 + 2a 3 ) /3 3 (a 3 + a 4 ) 
-4a 3 — a 4 — b H 5 1 5 



(A.9) 



and for k = 0, ±1. This solution was found in |8 
Its stability has been studied in 



(with the notation c 3 = a 3 , C4 = —0:4). 



21] . Although its physical metric is well behaved, the 



fiducial metric is degenerate, since both T and U do not depend on t. Moreover, T becomes 
imaginary for k = 1. 

Case lb: P2{u) = 0. This gives the same equation for u as in ( 15. bp . Therefore, one has 
b = u given by (ISTTj) . Eqs. flA.lblUA.lcp then yield 



T' = «aJ/«-l, 



T 



(A.10) 



and the integrability condition d\T = d^ t T reduces to (^^/ f k — 1J = \/ — 1 , which can 
be fulfilled only for k = — 1, when = sinh(r). As a result, the Stiickelberg scalars are 



T = u a(t) cosh(r), t/ = it a(t) sinh(r) 



(A.ll) 



The scale factor fulfills ( 1A.4I) with k = — 1 and with m 2 Tp = A, where A is the same as 
in (I5.13p . This solution was found in 10]. The fiducial metrics is non-degenerate, but the 



solution exists only for k = — 1. Moreover, it shows a non-linear instability 



io|. 



Let us now consider the case II, when the second factor in (1A.5I) vanishes, and so 

(1 + o?4 + 2a 3 )n — 3 — «4 — 3a 3 



(a 3 + a^ju — 1 — «4 — 2a 3 
The conservation condition (1A.3J) then reduces to 



a[u). 



(A.12) 



(b -a) {u-u x )P 1 {u) = 



(A.13) 



where 



U\ = —3 



a P 2 (u) 
a Pxiu)' 



P\{u) = 2(a 3 + a 4 )u — 4a 3 — 2a 4 — 2 



dP 2 {u) 
du 



(A.14) 
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Therefore, there are two options, either a = b or u = U\ (if P\(u) = 0, then a diverges). 
Case Ila: b = a. In this case Eqs. flA.ip yield T(t,r) = U(t,r), so that that the metric 
f^ v is degenerate. This also implies that a = b = 0, in which case one obtains from ( 1A.12I) 
u = (3«3 + «4 + 3)/(l + 2«3 + 04). The scale factor fulfills Eq. ( 1A.4[) with 



T ° = a 3 + 2 "3-«4 + 3 _ 

2a 3 + a 4 + l 
This solution has not been reported in the literature. 

Case lib: u — U\. Integrating this equation gives the relation between u and a, 

P 2 (n) = ^/a 3 , (A. 16) 



where A is the integration constant, so that u = u(t) and therefore a = a(t). Eqs. flA.laP 
and ( lA.lcP then yield 



T = \] a 2 + f3 2 f£, T'= jUMk : , (A.17) 

'o? + fi 2 fl 



with a(t) = aa and /3(t) = «a. The integrability condition d 2 r T = d 2 t T reduce to = Ba , 
where B is an integration constant. Using definitions of this transforms to 

B*a(u) Pl {u) 

uPxiu) - 3P 2 (u) ' 1 ; 

while Eq. flATbl) yields 

b= , Uf * . (A.19) 

VTTWf 2 

Since b depends both on t and r, there are three options to consider. The first two correspond 
to choosing either k = B = 0, or k = —1, B 2 = 1, in which cases b = u(t). The last option 
is to let b depend on r, but to set P2(u) = 0, in which case the coefficient in front of b in Tq 
vanishes. In each of these three cases, using (I A. 18[) to eliminate a in the Einstein equation 



( 1A.4|) and taking into account ( 1A. 16j) . one obtains an algebraic relation containing a and the 



matter density p. As a result, the solution exists only for fine-tuned p, which is unlikely to 
be physically interesting, so that we do not analyze this case any further. 
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